function [ratio_cumu,ratio_cumu_bangkok,ratio_cumu_northern,ratio_cumu_southern,ratio_cumu_central,ratio_cumu_northeastern]=...
    correct_ratio(comm_loc_model,comm_loc)
    
    load('region.mat');
    nyear=size(comm_loc,2);

    ratio_cumu(1:nyear-1,1)=0;
    ratio_cumu_bangkok(1:nyear-1,1)=0;
    ratio_cumu_northern(1:nyear-1,1)=0;
    ratio_cumu_southern(1:nyear-1,1)=0;
    ratio_cumu_central(1:nyear-1,1)=0;
    ratio_cumu_northeastern(1:nyear-1,1)=0;
    for i=1:nyear-1
        new_branch_cumu_model=find(comm_loc_model(:,1)==0 & comm_loc_model(:,i+1)==1);
        new_branch_cumu_data=find(comm_loc(:,1)==0 & comm_loc(:,i+1)==1);
        ratio_cumu(i)=length(intersect(new_branch_cumu_model,new_branch_cumu_data))/max(length(new_branch_cumu_model),length(new_branch_cumu_data));
    
    
        new_branch_cumu_model_bangkok=intersect(index_bangkok,find(comm_loc_model(:,1)==0 & comm_loc_model(:,i+1)==1));
        new_branch_cumu_data_bangkok=intersect(index_bangkok,find(comm_loc(:,1)==0 & comm_loc(:,i+1)==1));
        ratio_cumu_bangkok(i)=length(intersect(new_branch_cumu_model_bangkok,new_branch_cumu_data_bangkok))/max(length(new_branch_cumu_model_bangkok),length(new_branch_cumu_data_bangkok));
        
        new_branch_cumu_model_southern=intersect(index_southern,find(comm_loc_model(:,1)==0 & comm_loc_model(:,i+1)==1));
        new_branch_cumu_data_southern=intersect(index_southern,find(comm_loc(:,1)==0 & comm_loc(:,i+1)==1));
        ratio_cumu_southern(i)=length(intersect(new_branch_cumu_model_southern,new_branch_cumu_data_southern))/max(length(new_branch_cumu_model_southern),length(new_branch_cumu_data_southern));
    
        new_branch_cumu_model_northern=intersect(index_northern,find(comm_loc_model(:,1)==0 & comm_loc_model(:,i+1)==1));
        new_branch_cumu_data_northern=intersect(index_northern,find(comm_loc(:,1)==0 & comm_loc(:,i+1)==1));
        ratio_cumu_northern(i)=length(intersect(new_branch_cumu_model_northern,new_branch_cumu_data_northern))/max(length(new_branch_cumu_model_northern),length(new_branch_cumu_data_northern));
        
        new_branch_cumu_model_central=intersect(index_central,find(comm_loc_model(:,1)==0 & comm_loc_model(:,i+1)==1));
        new_branch_cumu_data_central=intersect(index_central,find(comm_loc(:,1)==0 & comm_loc(:,i+1)==1));
        ratio_cumu_central(i)=length(intersect(new_branch_cumu_model_central,new_branch_cumu_data_central))/max(length(new_branch_cumu_model_central),length(new_branch_cumu_data_central));

        new_branch_cumu_model_northeastern=intersect(index_northeastern,find(comm_loc_model(:,1)==0 & comm_loc_model(:,i+1)==1));
        new_branch_cumu_data_northeastern=intersect(index_northeastern,find(comm_loc(:,1)==0 & comm_loc(:,i+1)==1));
        ratio_cumu_northeastern(i)=length(intersect(new_branch_cumu_model_northeastern,new_branch_cumu_data_northeastern))/max(length(new_branch_cumu_model_northeastern),length(new_branch_cumu_data_northeastern));
    end



end